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Ideally, strong non-linearities could be used to implement quantum gates for photonic qubits by 
well controlled two photon interactions. However, the dependence of the non-linear interaction on 
frequency and time makes it difficult to preserve a coherent pulse shape that could justify a single 
mode model for the time-frequency degree of freedom of the photons. In this paper, we analyze the 
problem of temporal multi-mode effects by considering the pulse shape of the average output field 
obtained from a coherent input pulse. It is shown that a significant part of the two photon state 
transformation can be derived from this semi- classical description of the optical non-linearity. The 
effect of a non-linear system on a two photon state can then be determined from the density matrix 
' dynamics of the coherently driven system using input-output theory. As an example, the resonant 

(^*) , non-linearity of a single two level atom is characterized. The results indicate that the most efficient 

■ non-linear effect may not be the widely studied single mode phase shift, but the transfer of one of 

' the photons to an orthogonal mode distinguished by its temporal and spectral properties. 

PACS numbers: 42.50.Ex, 42.65.Sf, 42.50.Ar, 42.50.Ct 
^ ■ I. INTRODUCTION 



•j^ Photonic qubits based on the polarization or the transverse mode structure of single photons are an attractive 

, candidate for the implementation of quantum information protocols because it is relatively easy to establish and 
^ ' maintain single qubit coherences by conventional linear optics. Hovifever, the realization of well-controlled two photon 
interactions remains a challenging problem on the road to larger networks and more efficient operations. Early on, it 
has been suggested that non-linear materials might provide the interaction necessary to couple pairs of photons , 
I and sufficiently strong non-linear effects were demonstrated experimentally using single atoms in a high-finesse cavity 
• Q . Recent advances in solid state cavity designs seem to be putting the prospect of integrated devices implementing 
\^ I optical non-hnear quantum gates within the reach of present technological capabilities 0, SSBQ- However, some 
C — ■ . quite fundamental problems still need to be addressed before the proper functions of a quantum gate can be realized. 

' In particular, it is necessary to preserve single mode coherence, not only in polarization and in transverse mode 
' I structure, but also in the time-frequency domain. The latter problem is fundamentally linked to the dynamics of 
optical non-linearities in time and space [1, [l^, [Tl| ^ 1 ^ 13 Jl^. [Tsl . [T^ . \v^ . and has recently been identified as a 
critical problem in the realization of quantum gates IsT 13^" 2g| . Initial attempts at optimizing the pulse shape and 
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the pulse durations focused on the possibility of obtaining a large phase flip, represented by a negative two photon 
I amplitude in a single intended target mode [Til, [13. In that context, Koshino and Ishihara pointed out that the two 
I> ■ photon amplitude can be evaluated from the semi-classical response of the system [ll|. This result suggests that a 
more detailed analysis of the relation between semi-classical field expectation values and the transformation of two 
^ j , photon wavefunctions may be possible. Such an analysis could provide both a more thorough foundation for the 
■ evaluation of experimental results like the ones recently reported in [13, [l^ and a more detailed characterization 
and classification of the spectral and temporal features observed in highly non-linear devices. In the following, we 
therefore present a systematic analysis of the relation between the average field response obtained from a coherent 
input and the transformation of the two photon component of the quantum mechanical wavefunction. 

In section |lTl a quantum mechanical formulation of the semi-classical non-linear response is developed using field 
operators. In section Hill the non-linear interaction is expanded in terms of its effects on photon number states of 
suitably defined modes. It is then possible to express the semi-classical output in terms of the wavefunctions of 
these modes. In section IIV[ it is shown how the matrix elements of the photon number expansion can be derived 
from the overlap integrals of the semi-classical result. According to the results of sections [Til to [TVl the non-linear 
transformation has two distinct effects: a change in phase and amplitude of the two photon component of the linear 
output pulse and a change in pulse shape represented by the transfer of a single photon to an orthogonal mode. In 
section I VII the magnitude of these effects is investigated for the case of a resonant two-level system. It is shown 
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that the conditional transfer of one photon to an orthogonal mode is much stronger than the non-linear phase shift, 
suggesting that it may be more efficient to use this effect in optical quantum gates. In section IVlIl the pulse shapes 
of the modes involved in the non-linear photon tranfer are presented and the dynamics of the effect is discussed. The 
conclusions are summarized in section rVIIIl 



II. SEMI-CLASSICAL CHARACTERIZATION OF A NON-LINEAR QUANTUM SYSTEM 

In order to analyze the quantum level effects of an optical non-linearity, we start by considering the effects of a 
quantum system on the pulse shape of a coherent input pulse. Such an input pulse can be described by a classical 
time dependent field amplitude a foin(^)j where b[jj^{t) is the wave function of a single mode light field normalized to 
one photon. Hence, the input quantum state is effectively a single mode coherent state | a). This state interacts 
with the non-linear quantum system, temporarily exciting it from its initial ground state. After this interaction, the 
system returns to the ground state, re-emitting any photons it might have absorbed in the process of the interaction. 




FIG. 1; Illustration of the semi-classical characterization of an optically non-linear system using coherent input light and 
homodyne detection of the output pulse shape. 

If decoherence effects can be neglected, the total effect of the interaction on the quantum state of the input pulse 
can then be written as a unitary transformation U acting only on the state of the light field. In general, this unitary 
transformation acts on the full continuum of modes in free space. A complete characterization of the output state would 
therefore require multi-photon coherences between all frequencies or times observed in the output [9, 16, 18]. However, 
we now simplify this analysis by considering only the expectation values of the output fields {b{t)). Experimentally, 
this would correspond to the averages of a homodyne measurement, as used in the pioneering work of Turchette et 
al. 0]. Fig. [2 shows a schematic illustration of this semi-classical characterization of a quantum level non-linearity. 
Theoretically, such averages are easily obtained from the density matrix dynamics of the system using input-output 
theory, as will be explained in more detail in section IVl 

In terms of the unitary operation U describing the effects of the non-linear system on the quantum state of the 
input light, the average output field is given by the expectation value of the field operator b{t) in the time domain, 

"out 

(t) = {b{t)) = {a I U^b{t)U I a) (1) 

We can then proceed to analyze the specific form of the output pulse &out(0 by expanding the output in terms of the 
input amplitude a. The non-linear effect arising from photon-photon interactions is given by the third-order term 
in this expansion. For sufficiently small amplitudes a, the semi-classical effect of photon-photon interactions on the 
transformation of an input pulse a bj^{t) can therefore be expressed in terms of the third-order non- linear respone, 

bontit) = ab^'\t) + a\afb^'\t), (2) 

where b^^^{t) is the pulse shape of the linear output, and is the pulse shape of the third-order non-linearity. 

If there is neither decoherence nor photon loss in the system, b'-^^^t) describes a normalized single output mode 
that characterizes the transformation of the one photon wavcfunction from ip[-^{t) — b^y^{t) to ^pi{t) — b'-^^t). On 
the other hand, is neither normalized nor does it represent a wavcfunction orthogonal to V'i(i)- However, it is 

possible to interpret the pulse shape b^'^''{t) as a linear superposition of a component proportional to 8'^^ another 
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component proportional to a normalized wavefunction ■02 {t) that is orthogonal to (t) . Since the wavefunctions i/'i (t) 
and ip2 {t) are orthogonal, they represent two distinct modes in the time- frequency continuum. It is therefore possible 
to quantize the light field by assigning separate annihilation operators ai and 02 to these orthogonal modes. The 
output pulse shape 6out (^) can then be expressed in terms of the two orthonormal wavefunctions ipi (t) and ip2 (t) and 
the expectation values (di) and (0,2) of their complex amplitudes, 

bont{t) = {ai)Mt) + {a2)Mt)- (3) 

The expectation values can be determined by analyzing the pulse shape functions 6(^^(t) and b'^^^{t) in eq.([2]). Thus, 
the semi-classical representation of the non-linearity in terms of field expectation values can be interpreted within a 
fully quantum mechanical model focussing on the two modes defined by the pulse shapes observed in the output field. 



III. QUANTUM MECHANICS OF THE PHOTON-PHOTON INTERACTION 

In the previous section, we have shown that the third-order non-linear response to a coherent input pulse can be 
described in terms of two quantized modes defined by the linear and non-linear parts of the average output field 
^out(^)- We can now use this two mode representation to formulate the matrix elements of the unitary transformation 
U that describe the transitions between the photon number states of these modes. If photon losses can be neglected, 
the unitary transformation preserves the total photon number and we can look at each subspace of fixed total photon 
number separately. In particular, the vacuum state will not be changed by the interaction with the system, so the 
effect of U on the zero photon subspace is simply given by 

U I vac.) =1 vac). (4) 

For low intensity fields {a 1), the expectation values {b{t)) of the field amplitudes are given by the coherence 
between the single photon wavefunction at t and the vacuum. Therefore, the linear part of the semi-classical pulse 
shape transformation is equivalent to the transformation of the single photon wavefunction. Specifically, the single 
photon input wavefunction tp[j^{t) = b[^{t) is transformed into the linear output mode wavefunction 'ipi{t) — 
In terms of the annihilation operators of the input mode a^^ and the linear output mode di, the effect of the unitary 
transformation U in the single photon subspace can therefore be expressed as 



U 



(flin I vac.)) = a\ I vac). (5) 



The transformation in the two photon subspace is more complicated, because the photon-photon interaction gener- 
ally entangles the wavefunctions of the two photons [l^. However, we know from eq.© that a significant part of the 
two photon output wavefunction can be described in terms of the modes di and 0,2 ■ We can therefore expand the effect 
of U in the two photon subspace in terms of these output modes. Specifically, the contributions to the expectation 
values of the annihilation operators in eq. ([3]) originate from an inner product of a one photon wavefunction generated 
by annihilating a photon from the output two photon wavefunction and the output one photon state. As a result, 
the components of the two photon output contributing to the averages in eq. ^ must have at least one photon in the 
linear output mode di. The components of the two photon output that contribute to 6out(^) are therefore (i) a state 
where both photons are in mode di and (ii) a state where one photon is in mode di and one photon is in mode d2. In 
addition, there is a third component | rest) that does not contribute to 6out(0 because it does not have any photons 
in the linear output mode di. Thus, the expansion of the two photon output state can be written as 

■'■ .;it -t I \ \ _ ^lliAt 't I \ I r-.^^t At 



U |^y|«in«in I vac.) J = — ^djd} | vac.) -|- Ci2d{d2 | vac.) -I- Cr \ rest), (6) 

where Cn, C12, and Cr are the two photon amplitudes characterizing the photon-photon interaction described by U. 
Cii describes the amplitude of obtaining both photons in the same mode as the single photon output, C12 describes 
a non-linear transfer of one photon to a well defined orthogonal mode, and Cr describes the amplitude of processes 
where both photons are transfered to other modes with shapes that cannot be identified using only the semi-classical 
output average &out(^)- 

We can now solve eq.([T|) by applying the unitary transformation U to the single mode coherent input state | a). 
In order to describe the third-order non-linearity, it is convenient to expand the coherent state, neglecting all terms 
that only contribute terms of fourth or higher order in a to the final field expectation value in eq.((T]). The coherent 
state can then be approximated by 

I "> « f 1 - ^) I vac.) + a f 1 - H!") al \ vac.) + ^df^d}^ | vac). (7) 
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The application of [/ to this input state results in an output state of 



U\a)^ (^1- I vac.) +a(^l- ^ | vac.) + ^ I ^^c.) + C^dlal \ vac.) + Cr \ rest)^ . (8) 

From this output state, we can obtain the expectation values of ai and 0,2, 

(01) = (V^ut I ai I Vout) =a + (Cii - l)a|a|^ 

C 

(02) = (V'out I 02 1 V'out) = (9) 

These expectation values establish the connection between the few mode formulation of &out(^) in eq.© and the 
expansion up to third-order in a given by eq.®. Specifically, the expression for foout(0 obtained by inserting the 
results given by eq.® into eq.® reads 

|2 f ir< n\ , '^12, 



bont(t) = aA{t)+a\a\' ^{Cn - l)^i(t) + ^^sWj • (10) 

Based on this relation, it is possible to obtain the two photon amplitudes Ci that characterize the unitary transfor- 
mation U by decomposing the non- linear output pulse shape of eq.© into its and V'2(i) components. 



IV. DERIVATION OF TWO PHOTON AMPLITUDES d FROM SEMI-CLASSICAL PULSE SHAPES 

The two photon amplitudes Ci characterize the essential properties of the non-linear system for applications as an 
optical quantum gate. In particular, Cn describes any non-linear phase shift, with Cn = —1 corresponding to the 
ideal controlled phase- flip needed for the implementation of a quantum controUcd-NOT [l,[lll,[m. On the other hand, 
C12 describes a well-controlled transfer of one photon to a new mode. Since the output is fully quantum coherent, this 
process may also be a suitable candidate for quantum information processing. Finally, the coefficient Cr represents a 
component of unknown coherence that may be interpreted as a quantitative representation of the dispersion problem 
discussed in [isj . 

Eq. (jlO|) shows how the two photon amplitudes Cn, C12, and Cr can be determined from the semi-classical descrip- 
tion of the non-linearity in terms of b'^^^t) and 6^'^)(f). Specifically, the amplitude of the V'l (^"Component in b^^^{t) 
is equal to Cn — 1 and the amplitude of the -02 (^)-component is equal to Ci2/V^- Since "01 (i) is equal to the linear 
component of the semi-classical output, the ^i{t)-a,urp\itude can be determined from the overlap integral of 

b''^\t) and 6^^-'(t). The two photon amplitude Cn is therefore given by 

Cii = 1+ / b^^\t)*b'^^\t)dt. (11) 



As mentioned above, this parameter describes non-linear phase shifts that do not change the pulse shape of the photon 
wavepackets. This kind of single mode phase shift has been the focus of most of the previous work on non-linear 
optical quantum gates. In fact, Cn is equivalent to the parameter previously introduced by Koshino and Ishihara 
to evaluate the performance of a quantum non-linearity based on a semi-classical result [ll| . Our analysis completes 
this approach by taking into account the details of the pulse shape b'^^^ (t) that describes the semi-classical effects of 
the third-order non-linearity. In terms of the quantum mechanical description, this results in the introduction of the 
additional output mode ip2{t) and the associated two photon amplitude C12 describing the transfer of one photon 
to this new mode. The magnitude of this amplitude can be determined by taking the total intensity of 6^^^ (t) and 
subtracting the intensity accounted for by Cn, 



IC12P = 2(^1 |6(3)(t)|2di_|Cn-lp) 



= 2 / \b^^\t)\^dt- 



6(i)(i)*6(3)(t)dt 



(12) 



The phase of C12 depends on the definition of "02 (i)- It can therefore always be set to zero. Since the wavefunction 
-02 (i) of the target mode can also be obtained from the pulse shape 6'-^-'(i), the amplitude C12 also describes a fully 
coherent process that may be used to implement well controlled non-linear operations on optical quantum states. 
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All effects that cannot be described by the two photon amplitudes Cn and C12 are summarized by the amplitude 
Cr- Since the two photon output wavefunction given by eq.® is normalized, can be obtained from 

|aP = l-|Cnp-|Ci2p (13) 

Because this amplitude is associated with a two photon wavefunction that cannot be described within the two mode 
expansion of section |TT1 it has to be regarded as a source of decoherence in any straightforward implementation of a 
non-linear quantum gate. It thus provides a quantitative measure of the dispersion problems raised in ^18^, and it 
is an interesting question to what extend the amplitude Cr can be minimized while retaining the desired non-linear 
effects described by Cn and C12. 




FIG. 2: Illustration of the quantum limit of the integral J b^^\t)*b^^^ {t)dt = Cu — 1 describing the overlap between the linear 
and the non-linear part of the semi-classical output pulse. 

Before moving on to specific examples of quantum level non-linearities, it may be interesting to consider the natural 
limits imposed on semi-classical third-order non-linearities by the amplitudes Cu and C12 of the two photon output 
component. While classical physics imposes no such limits, quantization makes it impossible to have non-linear effects 
for light field intensities far below the single photon level. In the present theory, this limit is expressed quantitatively 
as restrictions on the magnitude of b^'^^t) corresponding to the requirement that |Ciip -I- |Ci2p < 1. Considering 
only Cii, the overlap between b^'^^t) and the normalized linear output is limited by 



< 1. (14) 



Figure [5] shows this limit of the overlap integral between the non- linear and the linear output wavefunctions as a circle 
in the complex plane. Note that the real part of the overlap is always negative, indicating that all third-order non- 
linearities reduce the output amplitude. Moreover, phase shifts of \(j)\ < tt/2 are associated with a minimal amplitude 
reduction of (1 — cos0). Likewise, a change in pulse shape described by the photon interaction amplitude C12 imposes 
a minimum on the negative value of the overlap between b'^^^t) and 6'^^^(t). As C12 increases, the radius of the circle 
in figure [5] is reduced to y^l — |Ci2p, resulting in a minimal amplitude reduction of 

Re (^J b^'Htrb<^^Ht)d?j < - (1 - V/I-IC12P) . (15) 

Thus, any non-linearity includes a reduction of the coherent output amplitude in the linear mode ipiit) — b^^'^{t). In 
particular, the maximal non-linear change in pulse shape (C12 = 1) requires an overlap of —1 between 6(^^(t) and 

6(3) (t). 



V. DYNAMICS OF THE NON-LINEAR SYSTEM 



As shown in the previous section, the two photon amplitudes Ci can be determined from the output pulse shapes 
b^^\t) and b^^\t) obtained from a semi-classical analysis of the non-linear optical response. It is therefore possible 
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to determine a significant part of the transformation acting on a two photon wavepacket by solving the non-linear 
dynamics of the system in response to a specific input pulse shape. In general, the dynamics of absorption and 
emission in an optical system excited by a coherent light field pulse can be represented by the dynamics of the density 
matrix, where the excitation is driven by a term linear in the coherent field amplitude a. Initially, the system is in its 
ground state p^^\ Weak excitations are described by a linear response p^^\ such that the density matrix dynamics 
are approximately given by p{t) = The field 

^out(^) — ab^^^{t) emitted by the system is then determined 
by the dipole or field expectation values of p*-^^ according to input-output theory |2J|. To describe the non-linear 
response, the density matrix dynamics can be expanded to include higher orders of a. Specifically, the excitation 
of the density matrix caused by an input pulse of intensity |ap can be described by a term |app^^^. However, the 
dipole or field expectations of this term are zero in all symmetric systems, since this contribution to p is invariant 
under a change of sign in a. Hence the lowest order non-linearity is obtained from the third-order term ajapp*^'^^ that 
describes the effects of the excitations on the response of the system. As the example given in the following will show, 
this kind of expansion can be done sequentially, resulting in a fairly simple and straightforward integration of a series 
of linear response equations. 




FIG. 3: Schematic representation of an atom-cavity system. If the reflectivity R of the back mirror is close to 1, the cavity 
mediated emission determines the coupling F with the fleld reflected by the front mirror. 

One of the most important systems considered for the realization of non-linear optical quantum gates is a single atom 
in a cavity. As Turchette et al. demonstrated experimentally ^ , such a system exhibits a strong non- linear phase shift 
in the weak coupling regime, where the cavity dynamics can be adiabatically eliminated and the system response is 
described by the Bloch equations of a two level atom. The problem of evaluating the strength of the non-linear response 
in this system has recently attracted a lot of interest, largely motivated by the improved experimental possibilities 
due to advances in cavity design d, 0, M, B S d, d, 113 ■ Our theory permits us to analyze the temporal and spectral 
properties of this non-linear response in terms of realistic pulse shapes without the complications introduced by a full 
quantum analysis of the entangled two photon wavefunction. It should be noted that we achieve this without any 
approximations, simply by omitting the parts of the two photon wavefunction that do not contribute to the average 
field 

^out(^)- In fact, the third order solutions of the semi-classical Bloch equations discussed in the following can 
also be obtained from the two photon wavefunctions determined in if one of the two photons is projected onto the 
single photon wavefunction. Thus the fully quantum mechanical analysis of the two photon response gives the same 
results for foout(^) the semi-classical analysis. It only provides additional information about the elusive component 
I rest) , which usually reqires an infinite number of additional modes for its precise characteization due to the spectral 
and temporal entanglement of the two photon wavefunction [lo|. For the coefficient Cn, the exact correspondence of 
the two photon response with the semi-classical result was first pointed out by Koshino and Ishihara in where 
they suggested the application of a coherent amplitude of a = 1/^/2 to obtain the maximal semi-classical response 
corresponding to Cn. In addition to this evaluation of Cn, the more detailed analysis developed here allows us to 
identify the photon transfer amplitude C12 that describes the non- linear change of pulse shape in terms of a conditional 
transition between the modes in the two photon output. As will be shown below, this effect is actually much stronger 
than the anticipated phase shift and may therefore play a significant role in possible realizations of non-linear optical 
quantum gates. Finally, we can also determine the precise pulse shapes ipiit) ^^'^ "02(0 that characterize the optical 
modes involved in the output of the quantum operation. 
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Fig. [3] shows a schematic iUustration of a single atom cavity system and its input/output characteristics. An input 
pulse is incident on the front mirror of the cavity. The cavity field adiabatically couples the input pulse to 

the atom with a coupling strength corresponding to a cavity enhanced spontaneous emission rate of 2r. The atom 
is excited and re-emits the absorbed energy by dipole emission through the cavity. Interference between this dipole 
emission and the reflected input pulse then results in the output pulse 6out(^)- In tti6 following, we assume that the 
losses of the system are negligible, so that all of the photons emitted by the atom will be found in 6out(^)- The dipole 
response of the atom can be described by the well known optical Bloch equations [1, H^l , 

-r(&-){t)-i2V2f abinit) {&z)it) (16) 
-2r(^(a,)(t) + i^ 

+*V2f - a* , (17) 

where (T_ is the operator describing the atomic dipole and is the operator describing the excitation of the atom. 
In general, the output field can be determined from the input field and the corresponding dipole response ((t_) of the 
atom according to input-output theory ^g, "2^ . For a lossless system, the corresponding relation reads 

bondt) = a6in(t) +i\/2f(a_)(t). (18) 

In general, the dipole response (cr_)(t) is a non- linear function of the input amplitude a. In order to determine the 
linear and third-order response in a, we can now apply the procedure described at the beginning of this section to 
the Bloch equations (fT6|) and (fT7| . 

Initially, the atom is in the ground state, so the zero order density matrix is described by {(Tz)'-^\t) = — ^- The 
linear response of the atom then determines the first order dipole term a((T_)*^^^ according to the linear relaxation 
dynamics given by using (cTz)^^\t) = — i in eq. (|16p . 

|(a_)(i)(t) = -r(a_)(i)(t) + zV2fV(t). (19) 

This equation can be solved for any input pulse shape by simply integrating the linear response. In principle, the 
second order of the density matrix is obtained by using the first order result ((T_)^^) as part of the excitation term in 
the relaxation dynamics described by eq. (fT7|) . However, a comparison of eq.(|16|l and cq. (fT7|) shows that this equation 
is always solved by the absolute square of the first order dipole term, (o'i^'') = |((T_)'^^'p. We can then determine 
((T_)(3)(t) from eqdZlby using (azY^Ht) instead of (ct^ )(")(*), 

|(a_)(3)(t) =-r(<T_)(3)(t)-^2V2f6i„(t) |(a_)«(t)p. (20) 

This equation has the same form as eq. (|19|) and can therefore be solved by the same kind of integration, corresponding 
to the linear response of the dipole to a modified input pulse of —2b[j^\{a-)^^^'^. Specifically, the third-order response 
describes the reduction of the linear response by the saturation of the gradually excited atoms, as indicated by the 
negative sign of the modified input pulse. 

Having obtained the linear and third-order responses of the atomic dipole, we can now express the output field up 
to third order in a using input-output theory as given by eg. ljlSp . The result reads 

6out(t) = «(v(<)+^^/2f(a_)(l)(^)) + a\a\^ (zV2f (a_)(3)(t)) (21) 

Comparison with eq.© shows how the solutions for the expectation values of the atomic dipole define the linear 
and non-linear output pulse shapes 6^^^(t) and 6^'^^(t). We can thus obtain the semi-classical characterization of the 
optical non-linearity necessary for the determination of the two photon amplitudes Ci and the wavefunctions tjjiit) 
and tjj2 (t) describing the modes used for the quantization by solving the density matrix dynamics of the system up to 
third-order in the coherent excitation amplitude a. 

VI. PULSE DURATION DEPENDENCE OF TWO PHOTON AMPLITUDES 

In general, the two photon amplitudes Ci depend on the specific shape of the input pulse defined by &in(i). In 
the following, we focus on resonant pulses, since such pulses seem to be the most promising candidates for strong 
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non-linear effects [8|. For a fixed pulse shape, the non-linearity then depends only on the ratio between pulse duration 
T and the relaxation time l/F of the atomic dipole. By applying our method of analysis, we can determine this 
dependence of non-linear effects on the scaled pulse duration FT. 

To cover a sufficiently wide range of possible pulse shapes while keeping the calculations relatively simple and 
efficient, we have chosen the four pulse shapes given in table U The most notable difference between the pulse shapes 
is that the change of the field is not continuous for the rectangular and the rising exponential pulse, while it is 
continuous for the symmetric exponential pulse and the Gaussian pulse. Moreover, only the rising exponential pulse 
is not symmetric around its peak. It should thus be possible to get insights into the effects of discontinuities and 
symmetry on the non-linear transformation of the pulses. 



TABLE I: Definition of input pulse shapes h-^-^(t). 



Type of pulse 


Wavefunction for pulse duration T 


Rectangular pulse 




' 1/Vt for - T < t < 
else 


Rising exponential pulse 




' y/2/fcxp{t/T) for t<0 
for t > 


Symmetric exponential pulse 


bin(*) = ^ 


/2/Texp{-2\t\/T) 


Gaussian pulse 


binW = V2/(V^r)exp(-2tVr2) 



For any given pulse shape and pulse duration, the non-linear output wavefunctions can be determined by solving 
eas. (|19ll2"T1) . Using these results, it is then possible to determine the two photon amplitudes Cn, C12, and Cr from eqs. 
(Ilim3p . The dependence of the results on pulse duration T for each of the four input pulse shapes is shown in fig. U) 
All four pulse shapes show a very similar pulse duration dependence, indicating that the two photon amplitudes Ci do 
not depend much on the specific pulse shape. Significant non-linear effects can typically be observed for pulse durations 
between about 0.1/F and 100/F. For shorter pulses (FT < 0.1), the bandwidth is too broad for resonant absorption, 
and for longer pulses (FT > 100), the photon density is too low for efficient non-linear interactions. Between these 
limits, the amplitude Cn drops from its linear value of 1 to negative values representing the non-linear phase flip 
originally proposed for use in non-linear optical quantum gates [2, S [lH. However, even the maximal values of the 
negative amplitude Cn represent only small fractions of the two photon output wavefunction. For the symmetric 
exponential input pulse (fig. [UJc)) and the Gaussian input pulse (fig. IHd)), the maximum is at about 0.2, while it 
is below 0.05 for the discontinuous pulse shapes of the rectangular input pulse (fig. [Ja)) and the rising exponential 
pulse (fig. m^b)). Thus, the non- linear phase flip is both limited in magnitude, and sensitive to discontinuities in the 
input pulse shape. 

On the other hand, the squared amplitude |Ci2p describing the probability of a non- linear transfer of exactly one 
photon to an orthogonal output mode V'2(i) has a peak value of about 2/3 for all four pulse shapes. This means that 
the conditional photon transfer is more efficient and less sensitive to pulse shape effects such as discontinuities. It may 
therefore be useful to consider this effect as an alternative option for the realization of non-linear optical quantum 
gates. 

VII. ANALYSIS OF OUTPUT PULSE SHAPES AT MAXIMAL PHOTON TRANSFER PROBABILITY 

Since the conditional photon transfer described by C12 is the strongest non-linear effect modifying the two photon 
output state, it may be useful to take a closer look at the coherent pulse shapes ipiit) ^-nd "02 (i) that determine 
about two thirds of the two photon output state. These pulse shapes can be obtained from the semi-classical results 
for 6(1) (t) and b'^^^t) using the amplitudes Cn and C12 and eqs. ([2]) and (jlOp. The maximal values of C12 are (a) 
|Ci2p = 0.66 at FT ~ 1.56 for rectangular input pulses, (b) |Ci2p = 2/3 at FT = 1 for rising exponential input 
pulses, (c) |Ci2p = 0.67 at FT = 0.78 for symmetric input pulses, and (d) |Ci2p = 0.64 at FT = 2 for Gaussian 
input pulses. For these pulse durations, the values of Cn are all vanishingly small, so the probability of finding both 
photons in the linear output mode is close to zero. 

Fig. [D shows the output pulse shapes at the maximal mode transfer amplitudes given above. Interestingly, the 
shapes of the non- linear target modes ip2{t) are somewhat similar in shape to the original input modes, while the 
wavefunctions ipi (t) are all delayed and change their sign to negative amplitudes after an initial positive peak (except 
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FIG. 4: Dependence of the squared two photon amplitudes \Ci\^ on pulse duration FT for (a) rectangular input pulses, (b) 
rising exponential input pulses, (c) symmetric exponential input pulses, and (d) Gaussian input pulses. Note that the pulse 
duration is given on a logarithmic scale. 



for the rising exponential input in fig. [n](b), where the amplitude is exactly zero until t = Q). This observation suggests 
a simple intuitive explanation for the non-linear mode transfer effect described by C12. If the pulse duration is perfectly 
matched to the absorption time, the atom is excited by one of the two photons for the entire pulse duration. Therefore, 
the second photon cannot be absorbed and passes the atom without the changes to its wavefunction otherwise induced 
by the linear dipole dynamics. If the transmitted mode "02 (i) is orthogonal to the linear response mode tpi{t), then 
the two photon output wavefunction will have exactly one photon in each of the two modes. A particularly interesting 
case may be that of the rising exponential input pulse shown in fig. [5] (b). Here, the input pulse and the output 
pulse are clearly separated in time. It may therefore be possible to distinguish the two modes by a sufficiently fast 
time-dependent gate. Specifically, any photon detected before t = indicates that a second photon was absorbed by 
the atom, since the linear output is exactly zero for t < 0. This effect could be used for a highly efficient elimination 
of multi-photon components in single mode quantum states. 



VIII. CONCLUSIONS 



We have presented an analysis of the relation between the average field amplitudes obtained from the density matrix 
dynamics of coherently driven systems and the transformation of two photon quantum states by a non-linear system. 
The results show how the spectral and temporal features of non-linear field transformations affect the performances 
of non-linear quantum gates operating on few photon states. It is thus possible to predict whether a given non-linear 
system can be used to implement a non-linear optical quantum gate operating on superpositions of quantum states 
with zero, one, and two photons, based on quantitative data of the intensity dependence of the coherent output field. 
In general, such data can be determined either experimentally, using coherent light inputs and homodyne detection, 
or theoretically, by solving the density matrix dynamics. 

The application to a resonant single-atom non-linearity shows that the most promising non-linearity may not be 
the widely investigated two photon phase shift, but a non- linear photon transfer process to a two photon output 
wavefunction where exactly one of the two photons is in a mode orthogonal to the single photon output mode. The 
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(a) (b) 





FIG. 5: Output pulse shapes tpiit) and tp2it) describing the linear output wavefunction and the target pulse of the conditional 
single photon transfer associated with C12. (a) shows the pulse shapes for a rectangular input pulse, (b) for a rising exponential 
input pulse, (c) for a symmetric exponential input pulses, and (d) for a Gaussian input pulse. 



investigation of the specific pulse shape shows that the effect may be understood in terms of the saturation of the 
atom when the linear output wavefunction ipi{t) is approximately orthogonal to the transmitted wavefunction ip2{t)- 
If the two wavefunctions can be separated by appropriate time-dependent gates, it may be possible to implement 
optical quantum gates based on this fundamental property of single atom non-linearities. 
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